Isospin-breaking corrections to superallowed Fermi /3-decay in isospin- and 
angular-momentum-projected nuclear Density Functional Theory 
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Background: The superallowed /3-decay rates provide stringent constraints on physics beyond the Standard 
Model of particle physics. To extract crucial information about the electroweak force, small isospin-breaking 
corrections to the Fermi matrix element of superallowed transitions must be applied. 

Purpose: Wc perform systematic calculations of isospin-breaking corrections to superallowed /3-decays and esti- 
mate theoretical uncertainties related to the basis truncation, time-odd polarization effects related to the intrinsic 
symmetry of the underlying Slater determinants, and to the functional parametrization. 

Methods: We use the self-consistent isospin- and angular-momentum-projected nuclear density functional theory 
employing two density functionals derived from the density independent Skyrme interaction. Pairing correlations 
are ignored. Our framework can simultaneously describe various effects that impact matrix elements of the Fermi 
decay: symmetry breaking, configuration mixing, and long-range Coulomb polarization. 

Results: The isospin-breaking corrections to the I = + ,T — 1 — > I = + ,T — 1 pure Fermi transitions are 
computed for nuclei from A=10 to A— 98 and, for the first time, to the Fermi branch of the J, T = 1/2 — ^ I, T = 1/2 
transitions in mirror nuclei from ^4=11 to A— 49. We carefully analyze various model assumptions impacting 
theoretical uncertainties of our calculations and provide theoretical error bars on our predictions. 

Conclusions: The overall agreement with empirical isospin-breaking corrections is very satisfactory. Using com- 
puted isospin-breaking corrections we show that the unitarity of the CKM matrix is satisfied with a precision 
better than 0.1%. 

PACS numbers: 21.10.Hw, 21.60.Jz, 21.30.Fe, 23.40.Hc, 24.80.+y 



I. INTRODUCTION 



By studying isotopes with enhanced sensitivity to fun- 
damental symmetries, nuclear physicists can test various 
aspects of the Standard Model in ways that are com- 
plementary to other sciences. For example, a possible 
explanation for the observed asymmetry between mat- 
ter and anti-matter in the universe could be studied by 
searching for a permanent electric dipole moment larger 
than Standard Model predictions in heavy radioactive 
nuclei that have permanent octupole shapes. Likewise, 
the superallowed /3-decays of a handful of rare isotopes 
with similar numbers of protons and neutrons, in which 
both the parent and daughter nuclear states have zero 
angular momentum and positive parity, are the unique 
laboratory to study the strength of the weak force. 

What makes these pure vector-current-mediated 
(Fermi) decays so useful for testing the Standard Model 
is the hypothesis of the conserved vector current (CVC), 
that is, independence of the vector current on the nuclear 
medium. The consequence of the CVC hypothesis is that 
the product of the statistical rate function / and partial 
half-life t for the superallowed I — + ,T = 1 ^ / = 
+ , T — 1 Fermi /3-decay should be nucleus independent 



and equal to: 



K 



const , 



(1) 



where K/{hcf = 27r 3 Mn2/(m c c 2 ) 5 = 8120.2787(11) x 
10~ 10 GeV~ 4 s is a universal constant; Gy stands for the 
vector coupling constant for semi-leptonic weak inter- 
action, and is the nuclear matrix element of the 
isospin rising or lowering operator T±. 

The relation (1) does not hold exactly and must be 
slightly amended by introducing a set of radiative cor- 
rections to the /i-values, and a correction to the nuclear 
matrix element due to isospin-symmetry breaking: 



\M^\ 2 =2(1-S C ), 



(2) 



see Refs. [1-4] and references cited therein. Since these 
corrections are small, of the order of a percent, they can 
be approximately factorized and arranged in the follow- 
ing way: 

K 



M = ft{l + S' B ){1 + S NS - S c ) 



(3) 



with the left-hand side being nucleus independent. In 
Eq. (3), = 2.361(38)% stands for the nucleus- 

independent part of the radiative correction [5], 5'^ 
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is a transition-dependent (Z-dependent) but nuclear- 
structure- independent part of the radiative correction [2, 
5] , and <5ns denotes the nuclear-structure-dependent part 
of the radiative correction [2, 6]. 

In spite of theoretical uncertainties in the evalua- 
tion of the radiative and isospin-symmetry-breaking cor- 
rections, the superallowed /3-decay is the most pre- 
cise source of experimental information for determin- 
ing the vector coupling constant Gy, and provides us 
with a stringent test of the CVC hypothesis. In turn, 
it is also the most precise source of the matrix ele- 
ment V u d = Gy/G^ of the Cabibbo-Kobayashi-Maskawa 
(CKM) three-generation quark mixing matrix [2, 7- 
9]. This is so because the leptonic coupling constant, 
Gf,/(hc) 3 = 1.16637(1) x 1(T 5 GcW 2 , is well known 
from the muon decay [9]. 

The advantage of the superallowed /3-decay strategy re- 
sults from the fact that, within the CVC hypothesis, Vud 
can be extracted by averaging over several transitions in 
different nuclei. For precise tests of the Standard Model, 
only these transitions that have ft- values known with a 
relative precision better than a fraction of a percent are 
acceptable. Currently, 13 "canonical" transitions spread- 
ing over a wide range of nuclei from A = 10 to A = 74 
meet this criterion (have /i-values measured with accu- 
racy of order of 0.3% or better) and are used to evaluate 
the values of Gy and Kid [2]. 

In this work we concentrate on the isospin-breaking 
(ISB) corrections 5q, which were already computed by 
various authors, using a diverse set of nuclear models 
[2, 10-17]. The standard in this field has been set by 
Towner and Hardy (HT) [2] who used the nuclear shell- 
model to account for the configuration mixing effect, 
and the mean-field (MF) approach to account for a ra- 
dial mismatch of proton and neutron single-particle (s.p.) 
wave functions caused by the Coulomb polarization. In 
this study, which constitutes an extension of our earlier 
work [15], we use the isospin- and angular- momentum- 
projected density functional theory (DFT) . This method 
can account, in a rigorous quantum-mechanical way, for 
spontaneous symmetry-breaking (SSB) effects, configu- 
ration mixing, and long-range Coulomb polarization ef- 
fects. 

Our paper is organized as follows. The model is de- 
scribed in Sec. II. The results of calculations for ISB cor- 
rections to the superallowed + — > + Fermi transitions 
are summarized in Sec. III. The ISB corrections to the 
Fermi matrix elements in mirror-symmetric T = 1/2 nu- 
clei are discussed in Sec. IV. Section V studies a particu- 
lar case of the Fermi decay of 32 CI. Finally, the summary 
and perspectives are given in Sec. VI. 



II. THE MODEL 

The success of the self-consistent DFT approach to 
mesoscopic systems [18] in general, and specifically to 
atomic nuclei [19-21], has its roots in the SSB mech- 



anism that incorporates essential short-range (pairing) 
and long-range (spatial) correlations within a single de- 
formed Slater determinant. The deformed states pro- 
vide a basis for the symmetry-projected DFT approaches, 
which aim at including beyond-mean-field correlations 
through the restoration of broken symmetries by means 
of projection techniques [22]. 

A. Isospin- and angular-momentum-projected DFT 
approach 

The building block of the isospin- and angular- 
momentum-projected DFT approach employed in this 
study is the self-consistent deformed MF state \<p) that vi- 
olates both the rotational and isospin symmetries. While 
the rotational invariance is of fundamental nature and is 
broken spontaneously, the isospin symmetry is violated 
both spontaneously and explicitly by the Coulomb inter- 
action between protons. The strategy is to restore the 
rotational invariance, remove the spurious isospin mix- 
ing caused by the isospin SSB effect, and retain only the 
physical isospin mixing due to the electrostatic interac- 
tion [23, 24]. This is achieved by a rediagonalization of 
the entire Hamiltonian, consisting the isospin-invariant 
kinetic energy and Skyrme force and the isospin-non- 
invariant Coulomb force, in a basis that conserves both 
angular momentum and isospin. 

To this end, we first find the self-consistent MF state 
\<p) and then build a normalized angular- momentum- 
and isospin-conserving basis \ip; IMK; TT Z ) by using the 
projection method: 

| W IMK; TT Z ) = 1 PI.tAkIv), (4) 

V M(p;IMK;TT z 

where P£ T ^ and K stand for the standard isospin 
and angular-momentum projection operators: 

2T -I- 1 F* 

P? Z ,T Z =^~ J <% zTz {p T )R{p T ) sin p T dp T , (5) 

Pii.K = I Dti K {Sl)R{Sl) dSl, (6) 

where, R{Pt) — e~ l ^ TTy is the rotation operator about 
the y-axis in the isospace, £0p T (/?t) is the Wigner func- 
tion, and T z — (N — Z)/2 is the third component of the 
total isospin T. As usual, R(Cl) = e -nj, e -H3jy e -i».h 
is the three-dimensional rotation operator in space, Q = 
(a, /3, 7) are the Euler angles, D I MK (Q) is the Wigner 
function, and M and K denote the angular-momentum 
components along the laboratory and intrinsic z-axis, re- 
spectively [22, 25]. Note that unpaired MF states \ip) 
conserve the third isospin component T z ; hence, the one- 
dimensional isospin projection suffices. 

The set of states (4) is, in general, overcomplete be- 
cause the K quantum number is not conserved. This 
difficulty is overcome by selecting first the subset of lin- 
early independent states known as collective space [22], 
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which is spanned, for each / and T, by the so-called nat- 
ural states \ip; IM; TT Z )W [26, 27]. The entire Hamilto- 
nian - including the ISB terms - is rediagonalized in the 
collective space, and the resulting eigenfunctions are: 



\n;p;IM;T z ) = Yl 

i,T>\T z \ 



a^\cp;IM;TT z )^, (7) 



One of the advantages of the projected DFT as com- 
pared to the shell- model-based approaches [2, 28] is that 
it allows for a rigorous quantum-mechanical evaluation 
of the Fermi matrix element using the bare isospin oper- 
ators: 



t ± 4£(^±^)^t 1±1 , 



(9) 



fc=i 



where the index n labels the eigenstates in ascending or- 
der according to their energies. The amplitudes <4™^ 
define the degree of isospin mixing through the so-called 
isospin- mixing coefficients (or isospin impurities), deter- 
mined for a given nth eigenstate as: 



where Ti±i denotes the rank-one covariant one-body 
spherical-tensor operators in the isospace, see the discus- 
sion in Ref. [29] . Indeed, noting that each mth eigenstate 
(7) can be uniquely decomposed in terms of the original 
basis states (4), 

|m; V ; IM; T z ) = ]T /^'"^^ PL, K \<p), (10) 

K.T 

with microscopically determined mixing coefficients 
f < Kj m ' I \ the expression for the Fermi matrix element 
where the sum of norms corresponds to the isospin T between the parent state |m; <p; IM; T z ) and daughter 



An;<p)\2 

l ilT I ' 



(8) 



dominating in the wave function \n; (p; IM; T z ). 



state |n; ip; IM; T z ± 1) can be written as: 



(m; p; IM; T z \f T \n; ft IM; r 2 ±l)=i^^ f^ mJ) V^' / ^^I^,t^i t iPC±i > t,±i^IV'> 



TT' KK' 

± lg7r 2 2^ 2^ JKT JK'T' 
TT' KK' 



dn d 1 k * k , (fi) (p\Pi >Ti f i Tl p£ ±hTz±1 \ip), (ii) 



where tilde indicates the Slater determinant rotated in space: \ip) = = R(il)\ip). The matrix element appearing 

on the right-hand side of Eq. (11) can be expressed through the transition densities that are basic building blocks of 
the mult i- reference DFT [24, 30-32]. Indeed, with the aid of the identity 



pT rp f>T' _ r TM r TK rp p T' 

r K,M ±x ^ r M' ,K> — ^T'M' \n 2—1 ^-'T'K-v\v 1 \vr K _ vK i 1 



v=-\ 



which results from the general transformation rule for spherical tensors under rotations or isorotations, 



the matrix element entering Eq. (11) can be expressed as: 

( <P I Pt z , T z T l =p 1 Pt z ± 1 , T z ± 1 1 V' ) = Ct ' T f z ± 1 1 =p 1 ^T^T Z - v 1 v ( f I ^1 v P~T Z - v , T z ± 1 1 $ ) • 



(12) 



(13) 



(14) 



For unpaired Slater determinants considered here, the double integral over the isospace Euler angles in Eq. (11) can 
be further reduced to a one-dimensional integral over the angle 0t using the identity 



(15) 

which is the one-dimensional version of the transformation rule (13) valid for rotations around the Oz axis in the 
isospace. The final expression for the matrix element in Eq. (14) reads: 

•> / - IT' + 1 

{<p\T lv PZ M _ ViTt±1 W) = — 



f d/3 T sml3 T d^ 1Tz±1 (p\f lu e' tl3Tf y\^) 
Jo 



{-iy 



2T' + 1 



d@ T sin. p T d%! z _ liTx±1 . 

r 



d 6 rp 1 _ v {n,h,r), 



(16) 



where Pi v (£l, At, r) is the isovector transition density, and the double-tilde sign indicates that the right Slater 
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determinant used to calculate this density is rotated both 
in space as well as in isospace: \ip) = R(/3t)R(^1)\i/j)- The 
symbol Af(Q,/3 T ) = {(p\R(P T )R{Q,)\ip) denotes the over- 
lap kernel. 

Since the natural states have good isospin, the states 
(7) are free from spurious isospin mixing. Moreover, since 
the isospin projection is applied to self-consistent MF 
solutions, our model accounts for a subtle balance be- 
tween the long-range Coulomb polarization, which tends 
to make proton and neutron wave functions different, and 
the short-range nuclear attraction, which acts in an op- 
posite way. The long-range polarization affects globally 
all s.p. wave functions. Direct inclusion of this effect in 
open-shell heavy nuclei is possible essentially only within 
the DFT, which is the only no-core microscopic frame- 
work that can be used there. 

Recent experimental data on the isospin impurity de- 
duced in 80 Zr from the giant dipole resonance 7-decay 
studies [33] agree well with the impurities calculated us- 
ing isospin-projected DFT based on modern Skyrme- 
force parametrizations [16, 23]. This further demon- 
strates that the isospin-projected DFT is capable of cap- 
turing the essential piece of physics associated with the 
isospin mixing. 

B. The choice of Skyrme interaction 

As discussed in Ref. [24], the isospin projection tech- 
nique outlined above does not yield singularities in en- 
ergy kernels; hence, it can be safely executed with all 
commonly used energy density functionals (EDFs) . How- 
ever, as demonstrated in Ref. [14], the isospin projection 
alone leads to unphysically large isospin mixing in odd- 
odd N = Z nuclei. It has thus been concluded that - 
in order to obtain reasonable results - isospin projection 
must be augmented by angular-momentum projection. 
This not only increases the numerical effort, but also 
brings back the singularities in the energy kernels [14] and 
thus prevents one from using the modern parametriza- 
tions of the Skyrme EDFs, which all contain density- 
dependent terms [34]. Therefore, the only option [14] is 
to use the Hamiltonian-driven EDFs. For the Skyrme- 
type functionals, this leaves us with one choice: the SV 
parametrization [35] . In order to better control the time- 
odd helds, the standard SV parametrization must be aug- 
mented by the tensor terms, which were neglected in the 
original work [35]. 

This density-independent parameterization of the 
Skyrme functional has the isoscalar effective mass as low 
as ~ 0.38, which is required to reproduce the actual 
nuclear saturation properties. The unusual saturation 
mechanism of SV has a dramatic impact on the over- 
all spectroscopic quality of this force, impairing such key 
properties like the symmetry energy [14], level density, 
and level ordering. These deficiencies also affect the cal- 
culated isospin mixing, which is a prerequisite for real- 
istic estimates of 5c ■ In particular, in the case of 80 Zr 



discussed above, SV yields «c w 2.8%, which is consider- 
ably smaller than the mean value of etc ~ 4.4 ± 0.3% ob- 
tained by averaging over nine popular Skyrme EDFs in- 
cluding the MSkl, SkO', SkP, SLy4, SLy5, SLy7, SkM*, 
SkXc, and SIII functionals, see Ref. [16] for further de- 
tails. Even though the ISB corrections 5c are primarily 
sensitive to differences between isospin mixing in isobaric 
analogue states, the lack of a reasonable Hamiltonian- 
based Skyrme EDF is probably the most critical defi- 
ciency of the current formalism. 

The aim of this study is (i) to provide the most reliable 
set of the ISB corrections that can be obtained within the 
current angular-momentum and isospin-projected single- 
reference DFT, and (ii) explore the sensitivity of re- 
sults to EDF parameters, choice of particle-hole configu- 
rations, and structure of time-odd fields that correlate 
valence neutron-proton pairs in odd-odd N = Z nu- 
clei. In particular, to quantify uncertainties related to 
the Skyrme coupling constants, we have developed a new 
density-independent variant of the Skyrme force dubbed 
hereafter SHZ2, see Table I. 

The force was optimized purposefully to properties of 
light magic nuclei below 100 Sn. The coupling constants 
to, ti,t2,XQ of SHZ2 were found by means of a \ 2 mm - 
imization to experimental [36] binding energies of five 
doubly-magic nuclei: 16 0, 40 Ca, 48 Ca, 56 Ni, and 100 Sn. 
The procedure reduced the x 2 from ~ 6.0 for SV set 
to ~ 3.6 for SHZ2. Most of the nuclear matter char- 
acteristics calculated for both sets are similar. It ap- 
pears, however, that the fit to light nuclei only weakly 
constrains the symmetry energy. The bulk symmetry en- 
ergy of SHZ2 is a sym w 42.2 MeV, i.e., it overestimates 
the accepted value a sym « 32 ± 2 MeV by almost 30%. 
While this property essentially precludes using SHZ2 in 
detailed nuclear structure studies, it also creates an in- 
teresting opportunity for investigating the quenching of 
ISB effects due to the large isospin-symmetry-restoring 
components of the force. 



TABLE I. Skyrme parameters U,Xi (i = 0,1,2,3), and W 
of SV [35] (second column) and SHZ2 (third column). The 
last column shows relative changes of parameters (in per- 
cent). Both parametrizations use the nucleon-mass parameter 
of ti 2 /2m = 20.73 MeV fm 2 . Parameters not listed are equal 



to zero. 


param. 


SV 


SHZ2 


change (%) 


to 


-1248.290 


-1244.98830 


-0.26 


ti 


970.560 


970.01156 


-0.06 


t 2 


107.220 


99.50197 


-7.20 


Xo 


-0.170 


0.01906 


-111.21 


W 


150 


150 
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C. Numerical details 

All calculations presented below were done by using 
the code HFODD [26, 37], version (2.48q) or higher, 
which includes both the angular-momentum and isospin 
projections. In order to obtain converged results for 
isospin mixing with respect to basis truncation, in our SV 
calculations we used N = 10 harmonic oscillator (HO) 
shells for A < 40 nuclei, 12 shells for 40 < A < 62 nuclei, 
and 14 shells for A > 62 nuclei. In SHZ2 test calcula- 
tions, we took N = 10 shells for A < 40 and N = 12 
shells for A > 40 nuclei. 

For the numerical integration over the Euler angles 
in space and isospace (a, j8, 7; 0t) we used the Gauss- 
Tchebyschev (over a and 7) and Gauss-Legendre (over j3 
and Pt) quadratures. We took n a = np — n 7 = 20 and 
nfj T = 8 (or 10) integration points. This choice guaran- 
tees that the calculated values of 6c are not affected by 
the numerical integration error. 



III. ISB CORRECTIONS TO THE 
SUPERALLOWED 0+ -> 0+ FERMI 
TRANSITIONS 



Mean field 



four-fold degeneracy 
of the sp levels 

aligned configuration 
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V®31 or V®K 

anti-aligned configuration 
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Isospin projection 



v « it T=0 



T=l 



V 8 It 
V® It 



T=0 



FIG. 1. Left: two possible g.s. configurations of an odd-odd 
N—Z nucleus, as described by the conventional deformed MF 
theory. These degenerate configurations are called aligned 
(upper) and anti-aligned (lower), depending on what levels 
are occupied by the valence particles. The right panel shows 
what happens when the isospin-symmetry is restored. The 
aligned configuration is isoscalar; hence, it is insensitive to the 
isospin projection. The anti-aligned configuration represents 
a mixture of T=0 and T=l states. The isospin projection 
removes the degeneracy by lowering the T—Q level. 



The + — > + Fermi /3-decay proceeds between the 
ground state (g.s.) of the even-even nucleus \I = 0, T ~ 
1,T Z = ±1) and its isospin- analogue partner in the 
N = Z odd-odd nucleus, \I = 0,T « 1,T Z = 0). The 
corresponding transition matrix element is: 

a4 ±} = (I = 0,Tt*l,T z =±l\f ± \I = 0,Tt*sl,T z = 0). 

(17) 

The g.s. state |J = 0,T » 1,T Z = ±1) in Eq. (17) is 
approximated by a projected state 

\I = 0,Tnl,T z = ±1) = ]T ^It.ii^o V), (18) 

T>1 

where is the g.s. of the even-even nucleus obtained 
in self-consistent MF calculations. The state \ijj) is un- 
ambiguously defined by filling in the pairwise doubly de- 
generate levels of protons and neutrons up to the Fermi 
level. The daughter state \I = 0,T w l,T z = 0) is ap- 
proximated by 

\I = 0,T« 1,T Z = 0) = ^oV^V), (19) 

T>0 

where the self-consistent Slater determinant \ip) = |z7(g)7r) 
(or \v ® 7r)) represents the so-called anti-aligned configu- 
ration, selected by placing the odd neutron and the odd 
proton in the lowest available time-reversed (or signature- 
reversed) s.p. orbits. The s.p. configuration |z><g>7r) man- 
ifestly breaks the isospin symmetry as schematically de- 
picted in Fig. 1. The isospin projection from \ip) as ex- 
pressed by Eq. (19) is essentially the only way to reach 
the |T w 1, 1 = 0) states in odd-odd N = Z nuclei. 



A. Shape-current orientation 




FIG. 2. (Color online) Schematic illustration of relative orien- 
tations of shapes and currents in the three anti-aligned states 
\<px) {X), I^y) {Y), and \cpz) (Z) discussed in the text. The 
long (Oz), intermediate (Ox), and short (Oy) principal axes 
of the nuclear mass distribution are indicated by thick arrows. 
The odd-neutron (j„) and odd-proton (J n ) angular momen- 
tum oriented along the Ox, Oy, or Oz axes is shown by thin 
arrows. Note that in each case the total angular-momentum 
alignment, j v + j^, is zero. 

At variance with the even-even parent nuclei, the anti- 
aligned configurations in odd-odd daughter nuclei are 
not uniquely defined. One of the reasons, which was 
not fully appreciated in our previous work [15], is re- 
lated to the relative orientation of the nuclear shapes 
and currents associated with the valence neutron-proton 
pairs. In all signature-symmetry-restricted calculations 
for triaxial nuclei, such as ours, there are three anti- 
aligned Slater determinants with the s.p. angular mo- 
menta (alignments) of the valence protons and neutrons 
pointing, respectively, along the Ox, Oy, or Oz axes of 
the intrinsic shape defined by means of the long (Oz), 
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intermediate (Ox), and short (Oy) principal axes of the 
nuclear mass distribution. These solutions, hereafter re- 
ferred to as \(fx), \ fy), and \ifiz), are schematically illus- 
trated in Fig. 2. Their properties can be summarized as 
follows: 

• The three solutions are not linearly independent. 
Their Hartree-Fock (HF) binding energies may typ- 
ically differ by a few hundred keV. The differences 
come almost entirely from the isovector correlations 
in the time-odd channel, as shown in the lower 
panel of Fig. 3 for a representative example of 34 C1. 
Let us stress that these poorly-known correlations 
may significantly impact the ISB corrections, as 
shown in the upper panel of Fig. 3. 

• The type of the isovector time-odd correlations cap- 
tured by the HF solutions depends on the relative 
orientation of the nucleonic currents with respect 
to the nuclear shapes. Solutions oriented perpen- 
dicular to the long axis, \tpx) and \(fy), are usually 
similar to one another (they yield identical corre- 
lations for axial systems) and differ from \ipz), ori- 
ented parallel to the long axis, which captures more 
correlations due to the current-current time-odd in- 
teractions. 

• The three \T = 1, 1 = + ) states projected from the 
\(fx), |</?y), and \ipz) Slater determinants differ in 
energy by only a few tens of keV, see the lower panel 
of Fig. 3. Hence, energy-wise, they represent the 
same physical solution, differing only slightly due 
to the polarization effects originating from different 
components of the time-odd isovector fields. How- 
ever, since these correlations are completely absent 
in the even-even parent nuclei, they strongly im- 
pact the calculated 5c- The largest differences in 
Sq have been obtained for A = 34 and A = 74 
systems, see Fig. 3 and Tables II and III. 

• Symmetry-unrestricted calculations always con- 
verge to the signature-symmetry-conserving solu- 
tion \ipz) which, rather surprisingly, appears to be 
energetically unfavored (except for 18 F). In spite of 
our persistent efforts, no self-consistent tilted-axis 
solutions have been found. 




X Y Z 



shape-current orientation 

FIG. 3. Upper panel: the ISB corrections Sc for the 
+ — > + supcrallowed /3-decays 34 Ar— > 34 C1 (open circles) 
and 34 C1— > 34 S (full circles) determined for the shape-current 
orientations X, Y, and Z depicted schematically in Fig. 2. 
Lower panel: differences between the energies of the X and Y 
configurations, and the Z configuration in 34 C1. Full triangles 
correspond to the total HF energies and open triangles corre- 
spond to contributions from the time-odd isovector channel. 
Full dots show the total energy differences obtained for the 
angular-momentum and isospin-projected states. 



the spherical uf 7 / 2 and 7r/ 7 / 2 sub-shells. These configu- 
rations can be labeled in terms of the quantum number 
K as \vK®irK) with K =1/2, 3/2, 5/2, and 7/2. 

In the extreme shell- model picture, each of these states 
contains all the T = 1 and 1=0, 2, 4, and 6 com- 
ponents. Within the projected DFT picture, owing to 
configuration-dependent polarizations in time-odd and 
time-even channels, the situation is more complicated be- 
cause the Slater determinants \vK ® ttK) corresponding 
to different if -values are no longer degenerate. Conse- 
quently, for each angular momentum /, one obtains four 
different linearly-dependent solutions. Calculations show 
that in all 7=0 and T « 1 states of interest, the isospin 
mixing ac is essentially independent of the choice of the 
initial Slater determinant. In contrast, the calculated 
ISB corrections 5c and energies depend on K, see Fig. 4. 



B. Nearly degenerate AT-orbitals 

Owing to an increased density of s.p. Nilsson levels 
in the vicinity of the Fermi surface for nearly spherical 
nuclei, there appears another type of ambiguity in choos- 
ing the Slater determinants representing the anti-aligned 
configurations. Within the set of nuclei studied in this 
work, this ambiguity manifests itself particularly strongly 
in 42 Sc, where we deal with four possible anti-aligned MF 
configurations built on the Nilsson orbits originating from 



C. Theoretical uncertainties and error analysis 

Based on the discussion presented in Sees. Ill A and 
IIIB, the recommended calculated values of 5c for the 
superallowed + — > + /3-decay are determined by aver- 
aging over three relative orientations of shapes and cur- 
rents. Only in the case of A = 42, we adopt for 5c an 
arithmetic mean over the four configurations associated 
with different if-orbitals. 

To minimize uncertainties in ac and 5c associated with 
the truncation of HO basis in HFODD, we used different 
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FIG. 5. (Color online) The convergence of the ISB correction 
8 C to the 46 V-> 46 Ti superallowed /3-decay versus the number 
of oscillator shells considered. 



FIG. 4. Top: the ISB corrections to the 42 Sc-s> 42 Ca super- 
allowed /3-transition, calculated using SV (circles) and SHZ2 
(dots) forces by projecting the \vK <g> ttK) configurations in 
42 Sc for K = 1/2, 3/2, 5/2, and 7/2. From top to bot- 
tom, horizontal lines mark (i) the average ISB correction 
using SV (thick solid line); (ii) the value of Ref. [2] (dot- 
ted line); and (iii) 5c of Ref. [12]. Shaded regions mark the 
related uncertainties. Bottom: projected energies of states 
\K;I — + ,T ~ 1) in 42 Sc obtained from the configurations 
\vK ® kK), relative to the projected energy of the K = 1/2 
state. 



HO spaces in different mass regions, cf. Sec. II C. With 
this choice, the resulting systematic errors due the ba- 
sis cut-off should not exceed ~ 10%. To illustrate the 
dependence of Sc on the number of HO shells, Fig. 5 
shows the case of the superallowed 46 V— j> 46 Ti transition 
obtained by projecting from the \(fz) solution in 46 V. In 
this case, the parent and daughter nuclei are axial, which 
allows us to reduce the angular-momentum projection to 
one-dimension and extend the basis size up to N — 20 
HO shells. 

With increasing N, Sc increases, and asymptotically 
it reaches the value of 0.8096(12)%. This limiting value 
is about 6.7% larger than the value of 0.7587% obtained 
for N — 12 shells, that is, for a basis used to compute 
the 42 < A < 54 cases. For 62 < A < 74 nuclei, which 
were all found to be triaxial, we have used N = 14 shells. 
The further increase of basis size is practically impossible. 
Nonetheless, as seen in Fig. 5, a rate of increase of Sc 
slows down exponentially with N, which supports our 
10% error estimate due to the basis truncation. 

The total error of the calculated value of Sc includes 
the standard deviation from the averaging, a n , and the 
assumed 10% uncertainty due to the basis size: A(#c) = 
V^n + (0-l<5c) 2 - The same prescription for A(<5c) was 
also used in the test calculations with SHZ2, even though 
a slightly smaller HO basis was employed in that case. 

For A — 38 nuclei, our model predicts the unusually 
large correction Sc ~ 10%. The origin of a very differ- 
ent isospin mixing obtained for odd-odd and even-even 
members of this isobaric triplet is not fully understood. 



Most likely, it is a consequence of the poor spectroscopic 
properties of SV. Indeed, as a result of an incorrect bal- 
ance between the spin-orbit and tensor terms in SV, the 
2s 1/2 subshell is shifted up in energy close to the Fermi 
surface. This state is more sensitive to time-odd polariza- 
tions than other s.p. states around 40 Ca core, see Table 
I in Ref. [38]. The calculated equilibrium deformations 
(/3 2 ,7) of the T z ± 1 and T z = A = 38 isobaric triplet 
are very similar, around (0.090, 60°). In the following, 
the 38 K— >- 38 Ar transition is excluded from the calculation 
of the Vud matrix element. 



D. The survey of ISB corrections in 10 < A < 74 
nuclei 

The results of our calculations are collected in Tables II 
and III, and in Fig. 6. In addition, Fig. 7 shows the 

differences, Sq V ^ — <5c HT \ between our results and those 
of Ref. [2] . In spite of clear differences between SV and 
HT, which can be seen for specific transitions including 
those for A = 10, 34, and 62, both calculations reveal a 
similar increase of Sc versus A, at variance with the RPA 
calculations of Ref. [12], which also yield systematically 
smaller values. 

The ISB corrections used for further calculations of 
V u d are collected in Table II. Let us recall that our pref- 
erence is to use the averaged corrections and that the 
38 K— > 38 Ar transition has been disregarded. All other 
ingredients needed to compute the J 7 ?;- values, including 
radiative corrections 5^ and (5ns , are taken from Ref. [2] , 
and the empirical ft- values are taken from Ref. [4]. For 
the sake of completeness, these empirical ft- values are 
also listed in Table II. 

In the error budget of the resulting ^-values listed in 
Table II, apart from errors in the ft values and radiative 
corrections, we also included the uncertainties estimated 
for the calculated values of Sc, see Sec. IIIC. To con- 
form with HT, the average value Ft — 3073.6(12)s was 
calculated by using the Gaussian-distribution-weighted 
formula. However, unlike HT, we do not apply any fur- 
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TABLE II. Results of calculations for the superallowed transitions measured experimentally. Shown are: the empirical ft-values 
[4]; SV values of 5c calculated by projecting from the \fx), \<Py), and \ipz) Slater determinants, see Sec. Ill A; recommended 



mean 5q corrections (see Sec. Ill C) and the corresponding Tt- values; empirical Sq xp; corrections calculated by using Eq. (23) 
contributions coming from the individual transitions to the \ 2 budget in the confidence-level test; mean 5^ SHZ2 - ) 
and the corresponding .Ft-values. 



corrections 



Parent 
nucleus 


/* 

(<0 


°c 

(%) 


*(Y) 

°c 
(%) 


r(Z) 

°c 
(%) 


4 SV) 
(%) 


Tt 

(s) 


^(exp) 

(%) 




r(SHZ2) 
(%) 


Tt 

00 


T z = -1 : 
10 C 
14 
22 Mg 
34 Ar 


3041.7(43) 
3042.3(11) 
3052.0(70) 
3052.7(82) 


0.559 
0.303 
0.243 
0.865 


0.559 
0.303 
0.243 
0.997 


0.823 
0.303 
0.417 
1.475 


0.65(14) 
0.303(30) 
0.301(87) 

1.11(29) 


3062.1(62) 
3072.3(21) 
3080.5(75) 
3056(12) 


0.37(15) 
0.36(06) 
0.62(23) 
0.63(27) 


3.7 
0.8 
1.9 
3.1 


0.462(65) 
0.480(48) 
0.342(49) 
1.08(42) 


3067.8(49) 
3066.9(24) 
3079.2(72) 
3057(15) 



T z =0 : 

26 Al 
34 C1 

42 Sc 

46y 

50 Mn 
54 Co 
62 Ga 
74 Rb 



3036.9(09) 
3049.4(11) 
3047.6(12) 
3049.5(08) 
3048.4(07) 
3050.8(10) 
3074.1(11) 
3084.9(77) 



0.308 0.308 0.494 
0.809 0.679 1.504 



0.486 
0.460 
0.622 
0.925 
2.054 



0.486 
0.460 
0.622 
0.840 
1.995 



0.759 
0.740 
0.671 
0.881 
1.273 



0.370(95) 
1.00(38) 
0.77(27) 
0.58(14) 
0.55(14) 
0.638(68) 
0.882(95) 
1.77(40) 

7i = 
IKdl = 



3070.5(31) 
3060(12) 
3069.2(85) 
3074.6(47) 
3074.1(47) 
3074.0(32) 
3090.0(42) 
3073(15) 

3073.6(12) 
0.97397(27) 
0.99935(67) 



0.37(04) 
0.65(05) 
0.72(06) 
0.71(06) 
0.67(07) 
0.75(08) 
1.51(09) 
1.86(27) 

x 2 = 

Xd = 



0.0 
48.4 
0.5 
4.5 
3.1 
2.0 
44.0 
0.1 

112.2 
10.2 



0.307(62) 
0.83(50) 
0.70(32) 

0.375(96) 
0.39(13) 
0.51(20) 
0.49(11) 
0.90(22) 

Tt = 
\V ud \ = 



3072.5(23) 
3065(15) 
3071(10) 
3080.9(35) 
3079.2(45) 
3078.0(66) 
3102.3(45) 
3101(11) 

3075.0(12) 
0.97374(27) 
0.99890(67) 



TABLE III. Similar as in Table II, except for the unmeasured 



transitions. 


Parent 
nucleus 


(%) 


°c 

(%) 


*(Z) 

°c 

(%) 


r(SV) 

°c 
(%) 


,(SHZ2) 

°c 
(%) 


T z = -\: 












18 Ne 

26 Si 

ao S 


2.031 
0.399 
1.731 


1.064 
0.399 
1.260 


1.142 
0.597 
1.272 


1.41(46) 
0.47(10) 
1.42(26) 


0.72(30) 
0.529(77) 
0.98(21) 


T z = : 












18p 


1.819 


0.956 


0.987 


1.25(42) 


0.42(24) 


22 Na 


0.255 


0.255 


0.535 


0.35(14) 


0.216(86) 


30p 


1.506 


0.974 


1.009 


1.16(27) 


0.60(20) 


66 As 
70 Br 


0.956 
1.654 


0.925 
1.479 


1.694 
1.429 


1.19(38) 
1.52(18) 


0.64(12) 
1.10(52) 



ther corrections to Ft. This leads to |Kd| = 0.97397(27), 
which agrees very well with both the HT result [2], 

\ V ud T) \ = 0.97418(26), and the central value obtained 
from the neutron decay \V^\ = 0.9746(19) [9]. A 
survey of the |Vud| values deduced by using different 
methods is given in Fig. 8. By combining the value of 
|Vud| calculated here with those of |y u s| = 0.2252(9) and 
jVubj = 0.00389(44) of the 2010 Particle Data Group [9], 
one obtains 

|Kd| 2 + |Ks| 2 + |Kb| 2 = 0.99935(67), (20) 

which implies that the unitarity of the first row of the 
CKM matrix is satisfied with a precision better than 
0.1%. A survey of the unitarity condition (20) is shown 



in Fig. 9. 



It is worth noting that by using 5c values correspond- 
ing to the fixed current-shape orientations (|<px), |</?y), 
or |cpz)) instead of their average, one still obtains com- 
patible results for |Vud| and unitarity condition (20), sec 
Figs. 8 and 9. Moreover, the value of |Vud| obtained by 
using SHZ2 is only «0.024% smaller than the SV result, 
see Table II. This is an intriguing result, which indicates 
that an increase of the bulk symmetry energy - that tends 
to restore the isospin symmetry - is partly compensated 
by other effects. The most likely origin of this compensa- 
tion mechanism is due to the time-odd spin-isospin mean 
fields, which are poorly constrained by the standard fit- 
ting protocols of Skyrme EDFs [41-43]. For instance, if 
one compares the Landau-Migdal parameters character- 
izing the spin-isospin time-odd channels [41-43] of SV 
(.go = 0.57, g' = 0.31, gi = 0.46, g[ = 0.46) and SHZ2 
(go = 0.27, g' = 0.30, g x = 0.47, g[ = 0.47) one notices 
that these two functionals differ by a factor of two in the 
scalar-isoscalar Landau-Migdal parameter g . 



To illustrate the compensation mechanism related to 
the bulk symmetry energy and go, in Fig. 10 we plot 5c 
for the 14 0^ 14 N^ 14 C super-allowed 0+ 0+ transi- 
tions as functions of the bulk symmetry parameter ci sym 
for a set of SV-based Skyrme forces with systematically 
varied xo parameter. At a functional level, xq affects only 
two Skyrme coupling constants (see e.g. Appendix A in 
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FIG. 6. (Color online) ISB corrections to the superallowed 
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and (b) Tz = — >■ Tz = 1 transitions. Our adopted values 
from Table II (circles with error bars) are compared with ISB 
corrections from Refs. [2] (dots, shaded band marks errors) 
and [12] (triangles). 
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FIG. 7. Differences between the ISB corrections to the twelve 
accurately measured superallowed + — > + /3-transitions (ex- 
cluding A— 38) calculated in this work with SV and those of 
HT [2]. Circles and dots mark the T z = -1 T z = and 
T z = — > T z = 1 decays, respectively. The errors, calculated 
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FIG. 8. The matrix element \V u d\ deduced from the superal- 
lowed + — > + /3-decay (dots) for different sets of the Sc cor- 
rections calculated in: (a) Ref. [2]; (b) Ref. [12] with NL3 and 
DD-ME2 Lagrangians; this work, using the averaged values of 
Sc with (c) SV and (d) SHZ2 functionals. Gray circles show 
ISB corrections with SV calculated at fixed shape-current ori- 
entations X, Y, and Z (from left to right). Triangles mark 
values obtained from the pion-decay [39] and neutron-decay 
[9] studies. The open circle shows |Kid| deduced from the 
/3-decays in the T = 1/2 mirror nuclei [40] 
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FIG. 9. Similar as in Fig. 8 except for the unitarity condition 
(20). 



ISB correction in Fig. 10 exhibits a minimum indicating 
the presence of the compensation effect. Similar effect 
was calculated for the A = 34 transitions. Hence, it is 
safe to state that our exploratory calculations are indica- 
tive of the interplay between the symmetry energy and 
time-odd fields. 



c '=M\ 



1 3 

24 



h 

24 d 



2 + X 3 J Po : 



2 ~ ^3 ] Pa ■ 



(21) 
(22) 



The coupling constant Cf influences the isovector part of 
the bulk symmetry energy [44] while Cq affects go. The 



E. Confidence level test 

In this section, we present results of the confidence- 
level (CL) test proposed in Ref. [4] . The CL test is based 
on the assumption that the CVC hypothesis is valid up 
to at least ±0.03%, which implies that a set of structure- 
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FIG. 10. ISB corrections for 14 0^ 14 N and 14 N^ 14 C super- 
allowed + — > + /3-decays calculated for a set of SV-based 
Skyrme forces with systematically varied xq parameter, which 
affects the bulk asymmetry energy coefficient a sym and spin- 
spin fields. The arrows indicate x$ values corresponding to 
SV and SHZ2. 



dependent corrections should produce statistically con- 
sistent set of Ft- values. Assuming the validity of the cal- 
culated corrections (5ns [6] , the empirical ISB corrections 
can be denned as: 



s-(oxp) 



1 



Ft 



3NS 



ft(i+s f R y 



(23) 



By the least-square minimization of the appropriate x 2 , 
and treating the value of Ft as a single adjustable pa- 
rameter, one can attempt to bring the set of empirical 
values 6q xp ^ as close as possible to the set of 5 C . 

The empirical ISB corrections deduced in this way are 
tabulated in Table II and illustrated in Fig. 11. Ta- 
ble II also lists individual contributions to the x 2 bud- 
get. The obtained x 2 P er degree of freedom (rid = 11) 
is x 2 l n d — 10.2. This number is twice as large as that 
quoted in our previous work [15], because of the large 
uncertainty of 5c for the 34 C1— >- 34 S transition. Other 
than that, both previous and present calculations have 
difficulty in reproducing the strong increase for A = 62. 
Our x 2 /rid is also higher than the perturbative- model 
values reported in Ref. [4] (x 2 / n d — 1-5), shell model 
with Woods-Saxon (SM-WS) radial wave functions (0.4) 
[2], shell model with Hartree-Fock (SM-HF) radial wave 
functions (2.0) [3, 45], Skyrme-Hartree-Fock with RPA 
(2.1) [11] , and relativistic Hartree-Fock plus RPA model 
(RHF-RPA) [12], which yields x 2 /n d = 1.7. 

It is worth noting that after disregarding the two 
transitions that strongly violate the CVC hypothesis, 
34 C1— > 34 S and 62 Ga— > 62 As that, and then performing a 
new CL test for the remaining ten transitions (rid = 9), 
the normalized \ 2 drops to 1.9. Within this restricted set 
of data, the calculated \V u d\ = 0.97420(28) and unitarity 
condition 0.99978(68) almost perfectly match the results 
of Ref. [2]. 
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FIG. 11. Top: differences between the calculated ISB cor- 
rections and empirical values resulting from the CL test of 
Ref. [4]. The shaded area of width ±0.2% is added in order to 
better visualize the differences. Bottom: contributions from 
individual transitions to the \ 2 budget. Note the particu- 
larly large contributions from the 34 C1-> 34 S and B2 Ga -^ 62 Zn 
transitions that deteriorate the CL test. See text for details. 



F. ISB corrections in 78 < A < 98 nuclei 

Our projected DFT approach can be used to predict 
isospin mixing in heavy nuclei. The calculated ISB cor- 
rections and Q-values in 78 < A < 98 nuclei are listed 
in Table IV. The values of 5c are also shown in Fig. 12. 
Note that the predicted ISB corrections are here consid- 
erably smaller than those in A=70 and A=7A nuclei, see 
Tables II and III. For the sake of comparison, Fig. 12 also 
shows predictions of Ref. [46] for the 82 Nb— > S2 Zr tran- 
sition using the VAMPIR approach with either charge- 
independent Bonn A potential or charge-dependent Bonn 
CD potential. Note that our prediction is only slightly 
below the Bonn A result and significantly lower than the 
Bonn CD value. For the sake of completeness, it should 
be mentioned that our Q^-value of 10.379 MeV for this 
transition agrees well with Qp = 10.496 MeV (Bonn A) 
and 10.291 MeV (Bonn CD) calculated within the VAM- 
PIR approach. 

Our calculated values of 5c are in heavy nuclei consid- 
erably smaller than those obtained from a perturbative 
expression [4, 10, 47]: 



0.002645 



A 2/3 



(n + l){n + 1 + 3/2) (%), (24) 



where n and I denote the number of radial nodes and an- 
gular momentum of the valence s.p. spherical wave func- 
tion, respectively. Indeed, assuming the valence lpi/2 
state in A = 78, Eq. (24) yields £ c =1.54%. In heavier 
nuclei, where the spherical valence state is 0gg/2, Eq. (24) 
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TABLE IV. Results of calculations for the superallowed /3-decays in 78 < A < 98 nuclei: the isospin impurities in the parent and 
daughter nuclei; 8q for different shape-current orientations; averaged (recommended) Sq; calculated equilibrium deformations 
p2 and 7; and Q/3-values calculated here and estimated from the extrapolated masses of Ref. [36]. 









(p) 
(%) 


(D) 
a C 
(%) 


°c 
(%) 


°c 

(%) 


S (z) 
°c 

(%) 


j(SV) 

(%) 


4 SV) 


7 (sv) 
(deg) 


Q< th) 

(MeV) 


(MeV) 


T z — 




T z = 1 






















78 y 




78 Sr 


2.765 


0.976 


1.20 


1.19 


1.20 


1.20(12) 


0.004 


60.0 


10.471 


10.650* 


82 Nb 




82 Zr 


3.099 


1.408 


0.70 


0.91 


0.70 


0.77(13) 


0.036 


60.0 


10.379 


11.220* 


86 Tc 


— > 


86 Mo 


3.337 


1.518 


0.89 


0.89 


1.08 


0.95(13) 


0.122 


0.0 


10.965 


11.350* 


90 Rh 




90 Ru 


3.525 


1.608 


0.99 


0.99 


1.09 


1.02(11) 


0.161 


0.0 


11.465 


12.090* 


94 Ag 




94 Pd 


3.674 


1.689 


0.86 


0.86 


1.17 


0.96(18) 


0.136 


0.0 


11.896 


13.050* 


98 In 


-» 


98 Cd 


3.805 


1.771 


0.89 


0.89 


1.36 


1.05(25) 


0.057 


0.0 


12.343 


13.730* 
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FIG. 12. The ISB corrections to the superallowed 0+ ->• + 
transitions in heavy nuclei calculated in the present work 
(full dots). Vertical bars mark the ISB corrections to the 
82 Nb— » 82 Zr transition calculated in Ref. [46] by using the 
VAMPIR formalism with the charge- independent Bonn A and 
charge-dependent Bonn CD interactions. 



gives 5c values that increase smoothly from 1.30% in 
A = 82 to 1.64% in A = 98. 



IV. ISB CORRECTIONS TO THE FERMI 
MATRIX ELEMENTS IN MIRROR-SYMMETRIC 
T = 1/2 NUCLEI 

Transitions between the isobaric analogue states in 
mirror nuclei \T = 1/2,7, T z = -l/2)-)-|T = 1/2, 7, T z = 
+1/2) offer an alternative way to extract the J-Y-values 
[48] and V u d [40, 49]. Those transitions are mixed Fermi 
and Gamow- Teller, meaning that they are mediated by 
both the vector and axial-vector currents. Hence, the 
extraction of V u d requires - in addition to lifetimes and 
Q-values - measuring another observable, such as the 
beta- neutrino correlation coefficient, beta-asymmetry, or 
neutrino-asymmetry parameter [50, 51]. Moreover, the 
method depends on the radiative and ISB corrections to 
both the Fermi and Gamow- Teller matrix elements. In 
spite of these difficulties, current precision of determina- 
tion of V u d using the mirror-decay approach is similar to 



that offered by neutron-decay experiments [9, 40, 49], see 
also Figs. 8 and 9. 

Within our projected-DFT model, we performed sys- 
tematic calculations of ISB corrections to the Fermi ma- 
trix elements, 6^, covering the mirror transitions in all 
11 < A < 49 nuclei. Calculations were based on the 
Slater determinants corresponding to the lowest-energy, 
unrestricted-symmetry HF solutions. If the unrestricted- 
symmetry calculations did not converge, the projection 
was applied to the constrained HF solutions with imposed 
signature symmetry. These two types of solutions differ, 
in particular, in relative shape-current orientation, which 
also varies with A depending on the s.p. orbit occupied by 
an unpaired nucleon. It should be underlined, however, 
that the HF solutions corresponding to the /3-decay part- 
ners were always characterized by the same orientation 
of the odd-particle alignment with respect to the body- 
fixed reference frame. All calculations discussed in this 
section were performed by using the full basis of N = 12 
HO shells and the SV force. 

The obtained values of the ISB corrections to the Fermi 
transitions, 



<5X 



\(T= 1 -,I,T Z 



-r 2 i Ti 2 2 n 



(25) 

are collected in Table V and illustrated in Fig. 13. Since 
the calculations were performed in a relatively large ba- 
sis, the basis-cut-off-related uncertainty in could be 
reduced to approximately 5%, cf. Sec. IIIC. Except 
for one case, theoretical spins and parities of decaying 
states were taken equal to those found in experiment: 
^fth) = / fc Xp )(g- s -)- Only for A = 31, no I = 1/2 com- 
ponent was found in the HF wave function, and thus the 
lowest solution corresponding to 7J h s =5/2+ was taken 
instead. It should be mentioned that, owing to the poor 
spectroscopic quality of SV, the projected states corre- 
sponding to 7^ xp ^(g.s.) are not always the lowest ones. 
This situation occurs for A = 19, 25, and 45, where the 
lowest states have 7^ = 5/2+, 1/2+, and 3/2~, and 

the corresponding 5^ values are 0.308 %, 0.419%, and 
0.636%, respectively. A relatively strong dependence of 
the calculated ISB corrections on spin is worth noting. 
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TABLE V. Results of calculations for \T = 1/2,7, T z = -l/2)-t\T = 1/2,7, T z = +1/2} /3-decays between mirror nuclei: 
theoretical spin and parity assignments; isospin-mixing coefficients in the parent and daughter nuclei; ISB corrections calcu- 
lated in this work (asterisks denote results obtained within unrestricted-symmetry calculations); ISB corrections of Ref. [48]; 
quadrupolc equilibrium deformation parameters in the parent nuclei; and theoretical and experimental Qp values. 
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FIG. 13. Full circles: calculated values of the ISB corrections 
to the Fermi transitions in T = 1/2 mirror nuclei. Open 
circles with errors: results calculated by Severijns et al. [48]. 



The calculations also indicate an appreciable impact of 
the signature-symmetry constraint on 5^, in particular, 
in the p/-shell nuclei with A — 45, 47, and 49. A simi- 
lar effect was calculated for the + — > + transitions, see 



(5c-values at fixed shape-current orientations in Tables II 
and III. 



V. THE ISB CORRECTION TO THE FERMI 
DECAY BRANCH IN 32 C1 

The Vud values extracted by using diverse techniques 
including + — > + nuclear decays, nuclear mirror decays, 
neutron decay, and pion decay are subject to both exper- 
imental and theoretical uncertainties. The latter pertain 
to calculations of radiative processes and - for nuclear 
methods - to the nuclear ISB effect. The uncertainties 
in radiative and ISB corrections affect the overall preci- 
sion of Kid at the level of a few parts per 10 4 each [1, 51]. 
It should be stressed, however, that the ISB contribution 
to the error bar of V u d was calculated only for a single 
theoretical model (SM-WS). Other microscopic models, 
including the SM-HF [3], RH-RPA [12], and projected 
DFT [15], yield 6c corrections that may differ substan- 
tially from those obtained in SM-WS calculations. 

Inclusion of the model dependence in the calculated 
uncertainties is expected to increase the uncertainty of 
Vud ■ According to Ref. [52] the increase can reach even an 
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order of magnitude. In our opinion, a reasonable assess- 
ment of systematic errors (due to the model dependence) 
cannot be done at present, as it requires the assumption 
that all the nuclear structure models considered are ei- 
ther equally reliable or their performance can be graded 
in an objective way. 

A good way to verify the reliability of various models is 
to compare their predictions with empirically determined 
5q- Recently, an anomalously large value of 5c ~ 5.3(9)% 
has been determined from a precision measurement of 
the 7 yields following the /3-decay of J = 1+,T = 1 state 
in 32 CI to its isobaric analogue state (Fermi branch) in 
32 S [53]. This value offers a stringent test on nuclear- 
structure models, because it is significantly larger than 
any value of 5c in the A — An + 2 nuclei. The physi- 
cal reason for this enhancement can be traced back to a 
mixing of two close-lying / = 1 + states seen in 32 S at 
the excitation energies of 7002 keV and 7190 keV, respec- 
tively [54]. The lower one is the isobaric analogue state 
having predominantly T = 1 component while the higher 
one is primarily of T = character. 

The experimental value 5c ~ 5.3(9)% is consistent 
with the SM-WS calculations: 5 C « 4.6(5)%. In our 
projected-DFT approach, we also see fingerprints of the 
strong enhancement in 5c value in 32 CI as compared to 
other A = An + 2 nuclei. Unfortunately, a static DFT ap- 
proach based on projecting from a single reference state 
is not sufficient to give a reliable prediction. This is be- 
cause, as sketched in Fig. 14, there exist ambiguities in 
selecting the HF reference state. In the extreme isoscalar 
s.p. scenario, by distributing four valence protons and 
neutrons over the Nilsson s.p. levels in an odd-odd nu- 
cleus, one can form two distinctively different s.p. config- 
urations, see Fig. 14. 



FIG. 14. Schematic illustration of several possible mean-field 
configurations in the odd-odd Z — N — 2 (left) and even- 
even N — Z (right) nuclei. The pairs of proton (neutron) 
s.p. levels, labeled as ir and n (y and v), are assumed to be 
degenerated due to the intrinsic signature symmetry. The 
orbits v and ir carry the signature quantum number r = —i 
(a = 1/2) while v and 7f have r = i (a — —1/2). 



The total signature of valence particles determines the 
total signature of the odd-odd nucleus and, in turn, an 
approximate angular-momentum distribution in its wave 
function [55]; the total additive signature aT(mod2) = 
0(1) corresponds then to even (odd) spins in the wave 
function [56]. It is immediately seen that the anti-aligned 
configuration shown in Fig. 14 has «t = 0; hence, in the 
first approximation, it can be disregarded. In this sense, 
the reference wave function in 32 C1 (or, in general, in 
any N — Z = ±2 odd-odd nucleus) corresponds to the 
uniquely defined aligned state. As seen in Fig. 14, this 
does not hold for 32 S (or, in general, for any N = Z 
even-even nucleus), where one must consider two possi- 
ble Slater determinants having qt = 1, obtained by a 
suitable proton or neutron particle-hole excitation. 

The above discussion indicates that, contrary to tran- 
sitions involving the odd-odd N = Z nuclei studied in 
Sec. Ill, those involving even-even N — Z nuclei cannot 
be directly treated within the present realization of the 
model. To this end, the model requires enhancements in- 
cluding the configuration mixing (multi- reference DFT). 
Nevertheless, we have carried out an exploratory study 
by independently calculating two ISB corrections for the 
two configurations discussed above. These calculations 
proceeded in the following way: 

• We select the appropriate reference con- 
figurations which, in the present case, 
are: u[A, 5, 3, 3]tt[5, 6, 3, 3] in 32 C1 and tpv 
u[5, 5, 3, 3]7r[4, 6, 3, 3] and <pu: v[A, 6, 3, 3]7r[5, 5, 3, 3] 
in 32 S. The labels denote the numbers of 
neutrons and protons occupying the lowest 
Nilsson levels in each parity-signature block 
(7r,r) = (+,+i),(+, counting 
from the bottom of the HF potential well, as 
defined in Ref. [57]. 

• We determine the lowest \I n = 1+, T w 1, T z = -1} 
and \ipi-,r = 1+,T w 1,T Z = 0) (i = I, II) states by 
projecting onto subspaces of good angular momen- 
tum and isospin, and performing the if-mixing and 
Coulomb rediagonalization as described in Sec. II. 

• Finally, we calculate matrix elements of the Fermi 
operator T± and extract 5c- 

The resulting ISB corrections are 4f :) = 2.40(24)% 

and 5(f u>> — 4.22(42)% for the ipi and ipu configurations, 
respectively. As before, we assumed a 10% error due 
to the basis size (N = 10 spherical HO shells). Pro- 
jections from the same configurations cranked in space 
to (J y ) = lh (see discussion in Ref. [27]) leaves ISB 

corrections almost unaffected: 5^^ = 2.41(24)% and 

5(f u>> — 4.30(43) %. A simple average value would read 
5c = 3.4(10)%, which is indeed strongly enhanced as 
compared to the A = An + 2 cases. The obtained cen- 
tral value is smaller than both the empirical value and 
the SM-WS result. It is worth noting, however, that 
within the stated errors our mean value 3.4(10)% agrees 







(odd-odd: Z-N=2] 


even-even: Z=N ] 


anti-aligned configuration 

V V JT Jt 

v(jtjt)®jt or v(jrjr)i»jt 

aligned configuration 

V V Jt Jt 
v(jtJt)®Jt or v(jtJt)®Jt 


aligned configuration I 

V V Jt Jt 
(vv)jt®Jt or (vv)jt®Jt 

aligned configuration II 

V V Jt Jt 
(jtJt)V0V or (jtJt)v®V 



14 



with the SM-WS value 4.6(05)%. Whether or not the 
configuration-mixing calculations would provide a signif- 
icant enhancement is an entirely open question. 

VI. SUMMARY AND PERSPECTIVES 

Within the recently-developed unpaired projected- 
DFT approach, we carried out systematic calculations 
of isospin mixing effects and ISB corrections to the su- 
perallowed + — > + Fermi decays in 10 < A < 74 nuclei 
and /3-transitions between the isobaric analogue states in 
mirror T = 1/2 nuclei with 11 < A < 49. Our predictions 
are compared with empirical values and with predictions 
of other theoretical approaches. Using isospin-breaking 
corrections computed in our model, we show that the 
unitarity of the CKM matrix is satisfied with a precision 
better than 0.1%. We also provide ISB corrections for 
heavier nuclei with 78 < A < 98 nuclei that can guide 
future experimental and theoretical studies. 

We carefully analyze various model assumptions im- 
pacting theoretical uncertainties of our calculations: ba- 
sis truncation, definition of the intrinsic state, and config- 
uration selection. To assess the robustness of our results 
with respect to the choice of interaction, we compared SV 
results with predictions of the new force SHZ2 that has 
been specifically developed for this purpose. The compar- 
ison of SV and SHZ2 results suggest that ISB corrections 
are sensitive to the interplay between the bulk symmetry 
energy and time-odd mean-fields. 



While the overall agreement with the empirical values 
offered by the projected-DFT approach is very encour- 
aging, and the results are fairly robust, there is a lot 
of room for systematic improvements. The main disad- 
vantages of our model in its present formulation include: 
(i) lack of pairing correlations; (ii) lack of ph interac- 
tion (or functional) of good spectroscopic quality; (iii) 
the use of a single HF reference state that cannot accom- 
modate configuration mixing effects; (iv) ambiguities in 
establishing the HF reference state in odd and odd-odd 
nuclei caused by different possible orientations of time- 
odd currents with respect to total density distribution. 
The work on various enhancements of our model, includ- 
ing the inclusion of T — and T — 1 pairing within the 
projected Hartree-Fock-Bogoliubov theory, better treat- 
ment of configuration mixing using the multi-reference 
DFT, and development of the spectroscopic-quality EDF 
used in projected calculations, is in progress. 
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